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We consider the continuous-time random walk of a particle in a two-dimensional self-affine 
quenched random potential of Hurst exponent H > 0. The corresponding master equation is studied 
via the strong disorder renormalization procedure introduced in Ref. [C. Monthus and T. Garel, J. 
Phys. A: Math. Theor. 41 (2008) 255002]. We present numerical results on the statistics of the 
equilibrium time t eq over the disordered samples of a given size L x L for 10 < L < 80. We find 
an 'Infinite disorder fixed point', where the equilibrium barrier Y eq = hit eq scales as Tgg — L u 
where it is a random variable of order O(l). This corresponds to a logarithmically-slow diffusion 
CN 1 \r(t) - r(0) ~ (lnt) 1/H for the position f(t) of the particle. 

f_L| 1 I. INTRODUCTION 



Random walks and diffusion processes have been the subject of constant interest in mathematics and in physics 
during the last century, for two main reasons (i) they play a central role in probability theory, and present a large 
number of very nice mathematical properties (ii) they naturally appear in a great variety of situations in physics 
and in biology. It is thus important to understand the effects of quenched disorder on random walks : are the usual 
properties of random walks stable with respect to the presence of some disorder or inhomogeneity ? If not, what are 
the new properties induced by disorder? Among the various types of random walks in random media that have been 
considered in the past (see the reviews [U-Il] and references therein), we wish to focus here on the case of random 
walks in a two-dimensional self-affine random potential U(r). In a continuous framework, this model can be defined 
via the Langevin equation for the position f of the particle 

^=-W(rH#) (1) 
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where ff(t) is a white noise 



<N ■ < m(tH(t') >= 2T $(t - t')hi ( 2 ) 

> ; 

that would generate a Brownian diffusion in the absence of the random potential [/, and where the quenched random 
potential U(r) is self-affine with some Hurst exponent H 

o ' 
d 



[U(r)-U{f ')f ~ \r-r'\ 2H (3) 



G\ . The case of dimension d = 1 and Hurst exponent H = 1/2 corresponds to the random-force Sinai model where the 
logarithmically-slow behavior \r(t) — r(0)| ~ (\nt) 2 has been obtained via various exact methods (see for instance the 
^ . review [3j] and references therein) . Since this logarithmic behavior replaces the usual power-law behaviour x ~ \/t of 
the pure Brownian motion, the effect of disorder is extremely strong. In higher dimension d > 1, the model is not 
exactly solvable, but from scaling arguments on barriers, one still expects the analogous logarithmic scaling 0, 0| 



\m-m\~(^t) 1/H (4) 

However, to the best of our knowledge, this behavior has not been much tested, except in the preliminary unpublished 
numerical results of Pettini shown on Fig. 4.9 of the review [3]. The aim of this paper is to study a continuous-time 
lattice version of this model in dimension d — 2, via the strong disorder renormalization procedure introduced in [f| 
that can be applied to any master equation in arbitrary dimension. 

The paper is organized as follows. In section HU we recall the Weierstrass-Mandelbrot function method to generate 
numerically two-dimensional self-affine random potentials. In section IIII1 we explain how to use for the present 
case the strong disorder renormalization method introduced in [f| . In section IIV1 we present our numerical results 
concerning the statistics of the equilibrium time t eq over the disordered samples of a given size L x L. Our conclusions 
are summarized in section fVl 
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II. METHOD FOR GENERATING A TWO-DIMENSIONAL SELF-AFFINE RANDOM POTENTIAL 

Among the various methods that have been proposed in the literature to generate random functions of a given 
Hurst exponent (see the reviews @ and a comparative study of their performances in Q), we have found numerically 
that the method giving the best results for the correlation of Eq. [3] is the so-called Weierstrass-Mandelbrot function 
method, that we recall in this section. 
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FIG. 1: Random self-amne random potential U(x,y) on a square of size 100 x 100, obtained via the method of Eq. [6] for the 
value H — 0.5 of the Hurst exponent : (a) Example of one realization of the random potential U(x,y) . (b) Log-log plot of the 

correlation function C(r) = [U(ri) — U(r2)] 2 as a function of the distance r = \ri — fi| after averaging over angles and over 
disorder realizations : the slope is here 2H — 1, as it should for H — 0.5 (see Eq. |3J. 



A. Reminder on the Weierstrass-Mandelbrot function in dimension d = 1 

In dimension d=l, the Weierstrass-Mandelbrot function is defined by [1, H[ 

cos(2iT(f> n ) - cos(2tt7"x + 2wcf> Tl ) 
U ( x ) = 1^ Z^h ( 5 ) 

n=nmi n 

where the phases <j) n are independent and uniform in [0, 1], and where n m i n — — oo and n max = +oo. The function 
U(x) is fractal with Hurst exponent H on all scales : the frequencies 7™ are in geometric progression, in contrast with 
a Fourier transform that would correspond to an arithmetic progression. In the limit 7 — > 1, the discrete spectrum 
become dense and the function U (x) converge towards the fractional Brownian motion of exponent H . We refer to 
Q for more details on its mathematical properties and now discuss how to use it for numerical simulations. 

If one wishes to generate the potential U(x) at N discrete points x = 1,2, ...N, one has to choose the three 
parameters (n m i n , n max , 7) in the following way : 

(i) the maximal Fourier frequency u ma x associated to the lattice spacing Ax — 1 is uj max = l/(Aai) = 1. A 
convenient choice is thus n max = 0, corresponding to the maximal frequency ^ nmax = 1 in the sum of Eq. [5] 

(ii) the minimal Fourier frequency Lu m i„ associated to the sample size N is ui m in = 1/iV. Since we do not wish any 
periodicity of order N in the potential, we have to choose n m in such that the minimal frequency 7"™"* in the sum of 
Eq. [5] satisfies 7™™" -C l/N. 

(iii) finally, the parameter 7 determines the discretization of the frequency spectrum : the frequency 7™ have to be 
sufficiently dense. 

We now turn to the generalization to higher dimension. 



3 



B. Generalization to dimension d — 2 

To generalise Eq. [5] to higher dimension d > 1 , the idea [ltj [H| is to keep the principle of a sum over plane waves 
of various vectors k, where the modulus |fc| varies in geometric progression 7™, and where the angular part of k is 
uniformly distributed to insure isotropy. In dimension d — 2, this corresponds to [ToL [ill] 

1 cos(27T0„ !m ) - cos(27T7 ri (xcos27ra„ im + ysin27ra„ im ) + 27r0„ !m ) 
tf(z,y) = ^ -= ^ — — (6) 



m * — ' T 



n-n m ; n v max m — 1 



where the phases a nim and <p n ,m are independent and uniformly distributed in [0,1]. The new parameter m max 
fixes the number of wave vectors k of a given modulus |fc|. We have checked that this generalization proposed in 
[Tol [Tlj gives satisfactory numerical realizations of self-affine random potential (whereas the alternative generalization 
proposed in [l2|, [H| that are based on cartesian coordinates presents anisotropy) . 

Here we wish to generate the potential U(x,y) at N 2 discrete points where x = 1,2, ...N and y = 1,2, .., iV. One 
has then to choose the four parameters {n m i n ,n max ,m max ,j) to obtain good results for the two-point function of 
Eq. [3] for all pairs of points of the samples. For squares samples of linear size 10 < L < 80, we have found that the 
following set of parameters give satisfactory realizations of the potential U(x,y) for Hurst exponents 0.3 < H < 0.8 
■ n ma x = 0, n m i n = —150, m max = 100, 7 = 1.2. We show on Fig. [1] (a) an example of realization of the random 
self-affine random potential U(x, y) on a square of size 100 x 100, for the value H = 0.5 of the Hurst exponent. The 
corresponding correlation function is shown on Fig. [T] (b) on a log-log plot. 



III. STRONG DISORDER RENORMALIZATION PROCEDURE 

Strong disorder renormalization (see [3] for a review) is a ver y sp ecific type of RG that has been first developed 
in the field of quantum spins : the RG rules of Ma and Dasgupta [15| have been put on a firm ground by D.S. Fisher 
who introduced the crucial idea of "infinite disorder" fixed point where the method becomes asymptotically exact, 
and who computed explicitly exact critical exponents and scaling functions for one-dimensional disordered quantum 
spin chains [16j . This method has thus generated a lot of activity for various disordered quantum models [l4 |. and 
has been then successfully applied to various classical disordered dynamical models, such as random walks in random 
media fl7l [l8T | . reaction-diffusion in a random medium (l9j . coarsening dynamics of classical spin chains [20I. trap 
models |2l|. random vibrational networks [22j], absorbing state phase transitions (23[, zero range processes [2J| and 
exclusion processes [251 ] . 

For random walks in random media, the procedure introduced in Refs (l7l [l8| or in the recent work [26| are specific 
to the dimension d = 1. Here in dimension d = 2, the appropriate framework is the 'strong disorder renormalization' 
(RG) procedure introduced Q that can be defined for any master equation. In this section, we recall its principles 
for the present problem of a particle in a two-dimensional potential. 



A. Master Equation 

The master equation describing the evolution of the probability P t (f) to be at position r at time t can be written 

as 

= E P * ( f ') W{r'^ r) - P t (r) W out (r) (7) 

r ' 

where W (r 1 — >■ r) represents the transition rate per unit time from position r ' to r , and 

W out {r) = Y,W{f^r') (8) 

r' 

represents the total exit rate out of position r. 

For the two-dimensional random walk in the random potential U(f) at temperature T, we have chosen to consider 
the Metropolis dynamics defined by the transition rates 

W{?^r ') = <W'> min(l, e -^ f ')-^»/ T ) (9) 
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The first factor S^p >> means that the two positions are neighbors on the two-dimensional lattice, and the last factor 
ensures the convergence towards thermal equilibrium at temperature T via the detailed balance property 

e ~U(r)/T w (f j?') = e -U{? ')/T w (f / ^ ^ (1Q) 

B. Strong disorder renormalization rules 

For dynamical models, the aim of any renormalization procedure is to integrate over 'fast' processes to obtain 
effective properties of 'slow' processes. The general idea of 'strong renormalization' for dynamical models consists in 
eliminating iteratively the 'fastest' process. The RG procedure introduced in [5[ can be summarized as follows : 

(1) find the position f* with the largest exit rate W* ut 

Kut = Wont (O = max r - [W ou t (f)\ (11) 

(2) find the neighbors (fi, fg, f n ) of position f* , i.e. the surviving positions that are related via positive rates 
W(f* — > fi) > and W(fi — > f*) > to the decimated position f*. For each neighbor position f with i G (1, .., n), 
update the transition rate to go to the position fj with j € (1, .., n) and j^i according to 

W new (n fj) = W(f -> fj) + W{f t -^r*)x ti> (fj) (12) 

where the first term represents the 'old' transition rate (possibly zero), and the second term represents the transition 
via the decimated position f* : the factor W(fi — > f*) takes into account the transition rate to f* and the term 

rr {rj) ~ w out {f*) {U > 

represents the probability to make a transition towards fj when in f*. The In rates W(f* fi) and Wjfj — > f*) then 
disappear with the decimated position f* . Note that the rule of Eq. [T2lhas been recently proposed in [27| to eliminate 
'fast states' from various dynamical problems with two very separated time scales. The physical interpretation of this 
rule is as follows : the time spent in the decimated position f* is neglected with respects to the other time scales 
remaining in the system. The validity of this approximation within the present renormalization procedure is discussed 
in detail in Q. 

(3) update the exit rates out of the neighbors of f*, with i = 1, .., n either with the definition 

wzr (ft) = E w ' new ( f i -> o (14) 

r 

or with the rule that can be deduced from Eq. [12] 

W™?{fi) = W out (n) - W{f t n W{ C? n) (15) 

*V out 

The physical meaning of this rule is the following. The exit rate out of the position fi decays because the previous 
transition towards f* can lead to an immediate return towards fi. After the decimation of the position f* , this 
process is not considered as an 'exit' process anymore, but as a residence process in the position rj. This point is 
very important to understand the meaning of the renormalization procedure : the remaining positions at a given 
renormalization scale are 'formally' microscopic positions of the initial master equation (Eq. [7]), but each of these 
remaining microscopic position actually represents some 'valley' in position space that takes into account all the 
previously decimated positions. 

(4) return to point (1). 

We refer to [5| for more detailed explanations. In practice, the renormalized rates W(f — ¥ f) can rapidly become 
very small as a consequence of the multiplicative structure of the renormalization rule of Eq[T21 This means that the 
appropriate variables are the logarithms of the transition rates, that we will call 'barriers' in the remaining of this 
paper. The barrier B(f — > f ') from f to f is defined by 

B(f ->f') = -\nW(f ->r ') (16) 

and similarly the exit barrier out of position f is defined by 

B out (f) = - In W out (f) (17) 

A very important advantage of this formulation in terms of the renormalized transition rates of the master equation 
is that the renormalized barriers take into account the true 'barriers' of the dynamics, whatever their origin which 
can be either energetic or entropic. 



5 



C. Numerical details 



We have applied numerically these renormalization rules for square samples of size L 2 with 10 < L < 80 with 
a statistics of 33. 10 5 > n s (L) > 280 disordered samples. We have studied six values of the Hurst exponent in the 
interval 0.3 < H < 0.8. 



IV. STATISTICS OF THE EQUILIBRIUM TIME OF FINITE SYSTEMS 
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FIG. 2: (Color on line) Statistics of the equilibrium time t eq over the disordered samples of sizes L 2 : the appropriate variable 
is r eq = \nt e q (Eq. (a) the disorder- averaged value Y eq (L) and the width A(L) shown here for H = 0.8 scale with the 
same exponent ip (see Eq. I20[) . (b) the exponent ip of the width A(L) coincides with the Hurst exponent H, as shown here for 
H = 0.3,0.4,0.5,0.6,0.7,0.8 

In a finite system, the master equation of Eq. [7] satisfying the detailed balance condition of Eq. [10] will converge 
exponentially towards the equilibrium Boltzmann distribution. The characteristic time of this exponential convergence 
is called the equilibrium time t eq 

P t (r)-P eq {r) cx (18) 

t— ►oo 

Within the strong disorder renormalization procedure described in the previous section, this equilibrium time t eq 
of a given disordered sample is determined by the renormalized exit barrier 

r e , = lni e? (19) 

corresponding to the last decimation process where the two largest metastable valleys merge into a surviving valley 
corresponding to thermal equilibrium of the whole sample. We find that the disorder-averaged value T eq (L) and the 
width A(L) involve the same barrier exponent ip 

TT q (L) cx ifi 

L—^oo 

A(L) cx (20) 



L- 



as shown on Fig. d](a) for the value H = 0.8 of the Hurst exponent. Moreover, this exponent ip is equal, as expected 
[1, 0] j to the Hurst exponent H of the random potential 

ip = H (21) 

We show on Fig. [2] (b) the log-log plot of the width A(L) for various values H = 0.3, 0.4, 0.5, 0.6, 0.7, 0.8 of the Hurst 
exponent. Our conclusion is thus that the strong disorder renormalization procedure confirms the activated nature of 
the dynamics and the logarithmic-slow diffusion of Eq. [J] in dimension d = 2. 
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FIG. 3: (Color on line) Statistics of the equilibrium time t eq over the disordered samples of sizes L 2 for the case H — 0.8 : (a) 
Probability distribution Qhi^eq = lni e? ) for L — 10,20,30,40,50 (b) same data after the rescaling of Eq. 1221 and in log-scale 
to see the tails. 



We show on Fig. [3] (a) the probability distribution <2z,(r e(? = \nt eq ) over the disordered samples of size L 2 for 
various sizes L, for the case H = 0.8. The convergence towards a fixed rescaled distribution 

Q^-^O^^fW) (22) 
is shown on Fig. [3] (b) in log-scale to see the tails. 



V. CONCLUSION 

In this paper, we have shown that the strong disorder renormalization rules for master equations introduced in 5] 
are appropriate to study random walks in two-dimensional self-affine random potentials of Hurst exponent H > : 
we have found an 'Infinite disorder fixed point', where the equilibrium time t eq to reach equilibrium for samples of 
size L x L scales as \\it eq — L H u where u is a random variable of order O(l). This activated scaling found for the 
dynamics indicates that the strong disorder renormalization procedure becomes asymptotically exact in the limit of 
large times and large sizes [HJ • Our results confirm that the logarithmic-slow diffusion of Eq. |4] exists not only in 
dimension d = 1 where exact results can be obtained for the Sinai model case H = 1/2, but also in higher dimension, 
as shown here for d = 2. These conclusions have been recently checked via independent methods based on the exact 
calculation of the biggest relaxation time 28] or of some first-passage time [2!| . 
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